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Abstract 

A new Higgs-like boson with mass around 126 GeV has recently been discovered at the 
LHC. The available data on this new particle is analyzed within the context of two-Higgs 
doublet models without tree- level flavour-changing neutral currents. Keeping the generic 
Yukawa structure of the Aligned Two-Higgs Doublet Model framework, we study the 
implications of the LHC data on the allowed scalar spectrum. We analyze both the CP- 
violating and CP-conserving cases, and a few particular limits with a reduced number of 
free parameters, such as the usual models based on discrete Z2 symmetries. 



1 Introduction 



The ATLAS and CMS collaborations have recently announced the discovery of a new neutral 
boson, with a measured mass of 125.2 ± 0.3 ± 0.6 GeV [l] and 125.8 ± 0.4 ± 0.4 GeV |], 
respectively. The LHC data is compatible with the expected production and decay of the 
Standard Model (SM) Higgs boson, the most significant decay modes being H — > 77 and 
H ZZ^ £ + £~. The excess of events observed by ATLAS (CMS) has a (local) statistical 
significance of 6.1 a (6.9 a). Although the spin of the new particle has not been measured yet, 
the observed diphoton decay channel shows clearly that it is a boson with J 7^ 1, making very 
plausible the scalar hypothesis. Preliminary analyses of H — > ZZ — > A£ [3j[4] and H — > 77 [5j[6] 
events suggest indeed the assignment J p = + , though more statistics is still needed to give a 
definite answer. 

Additional (but less significant) evidence has been reported by the CDF and D0 collabo- 
rations |7j, which observe an excess of events in the mass range between 120 and 135 GeV (the 
largest local significance is 3.3 a). The excess seems consistent with a SM Higgs produced in 
association with a W ± or Z boson and decaying to a bottom-antibottom quark pair. 

While more experimental analyses are needed to assess the actual nature of this boson, the 
present data give already very important clues, constraining its couplings in a quite significant 
way. The stringent exclusion limits set previously on a broad range of masses provide also 
complementary information which is very useful to establish allowed domains for alternative 
new-physics scenarios. A SM Higgs boson has been already excluded at 95% CL in the mass 



ranges 0-122.5 and 127-600 GeV 5 8 14 



The new boson appears to couple to the known gauge bosons (W ± , Z, 7, g) with the 



strength expected for the SM Higgs 15 -22 , although a slight excess of events in the 27 decay 



channel, compared with the SM expectation, is observed by ATLAS and CMS [l]|2j. Moreover, 
its fermionic couplings seem compatible with a linear dependence with the fermion mass, scaled 



by the electroweak scale v ~ 246 GeV 23 . Thus, it has the properties expected for a Higgs- 
like particle, related with the spontaneous breaking of the electroweak symmetry. An obvious 
question to address is whether it corresponds to the unique Higgs boson incorporated in the 
SM, or it is just the first signal of a much richer scalar sector. 

The simplest modification of the SM Higgs mechanism consists in incorporating additional 



2 



scalar doublets, respecting the custodial symmetry, which can easily satisfy the electroweak 
precision tests. This leads to a rich spectrum of neutral and charged scalars, providing a 
broad range of dynamical possibilities with very interesting phenomenological implications. The 
minimal extension of the scalar sector with only one additional doublet contains five physical 
scalars: two charged fields H ± and three neutral ones h, H and A; thus, there are three possible 
candidates for the recently discovered neutral boson. If the scalar potential preserves the CP 
symmetry, h and H are CP-even, while A is CP-odd; in this case there are no AW + W~ and 
AZZ couplings at tree level, which makes the A possibility quite unlikely. 

Generic multi-Higgs doublet models give rise to unwanted flavour-changing neutral cur- 
rent (FCNC) interactions through non-diagonal couplings of neutral scalars to fermions. The 
tree-level FCNCs can be eliminated requiring the alignment in flavour space of the Yukawa ma- 
trices coupling to a given right-handed fermion (24] . The Aligned Two-Higgs Doublet Model 
(A2HDM) [25] results in a very specific structure, with all fermion-scalar interactions being 
proportional to the corresponding fermion masses. This leads to a rich and viable phenomenol- 



ogy 24-29 with an interesting hierarchy of FCNC effects, suppressing them in light-quark 
systems while allowing potentially relevant signals in heavy-quark transitions. The A2HDM 
constitutes a very general framework which includes, for particular values of its parameters, 



all previously considered two-Higgs doublet models (2HDMs) without FCNCs 30 31 , and 
incorporates in addition new sources of CP violation. 

In the following, we will analyze the recent discovery of a Higgs-like object within the 
A2HDM. We will study the different possible interpretations of the new boson, the correspond- 
ing experimental constraints on its couplings, and the implications for the remaining scalar 



spectrum. Previous analyses 32 -42 have only considered more specific scenarios based on dis- 



crete Z2 symmetries G3], i.e., the so called 2HDMs of types I 44 45 , II 45 46 , X (leptophilic 



or lepton specific), Y (flipped) [47 50 and inert 51 . The more general A2HDM framework 
opens a wide range of additional possibilities, which we will try to characterize keeping in mind 
the high-statistics data samples that the LHC is expected to deliver in the future, at higher 
energies. Two very recent works have already employed the A2HDM, in the limit of CP con- 



servation, to analyze the Higgs data 52 53 . We will compare our results in that limit and will 
also explore the consequences of allowing CP-violating phases, either in the scalar potential 
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(mixing of the three neutral scalars) or in the Yukawa couplings. 

Our paper is organized as follows: In section |2j we describe the theoretical framework 
adopted in our analysis, indicating the relevant couplings of the A2HDM scalars. In section 3 
we define the Higgs signal strengths, which are used to make contact with the experimental 
measurements. Section [4] presents our results and shows the scalar parameter ranges needed 
to explain the present data. Our conclusions are given in section [5j The appendices include a 
compilation of useful formulae as well as the statistical treatment and data used in this work. 



2 The Aligned Two-Higgs-Doublet Model 



The 2HDM extends the SM with a second scalar doublet of hypercharge Y = |. The neutral 
components of the scalar doublets <f> a (x) (a = 1,2) acquire vacuum expectation values that 
are, in general, complex: (O|0^(x)|O) = (0, v a e ida ). Through an appropriate U(l)y trans- 
formation we can enforce Q\ = 0, since only the relative phase 9 = 9 2 — 6\ is observable. It is 
convenient to perform a global SU(2) transformation in the scalar space (<pi,4>2) and work in 
the so-called Higgs basis ($i, $2), where only one doublet acquires a vacuum expectation value: 



(1) 







cos (3 sin (3 






sin (3 — cos j3 



with tan/3 —v 2 /vi. In this basis, the two doublets are parametrized as 
$1 = 



G+ 



V2 



$5 



H+ 
k (Si + iSs) 



(2) 



where G ± and G° denote the Goldstone fields and (0|#+|0) = (0|G+|0) = (0|G°|0) = (0|^|0) = 
0. Thus, $1 plays the role of the SM scalar doublet with v = yjv\ + v\ ~ (y/2G F y 1/2 = 
246 GeV. 

The physical scalar spectrum contains five degrees of freedom: the two charged fields H ± (x) 
and three neutral scalars $(x) = {h(x),H(x),A(x}}, which are related with the Si fields 
through an orthogonal transformation <Pi(x) = TZijSj(x). The form of the TZ matrix is fixed 
by the scalar potential, which determines the neutral scalar mass matrix and the correspond- 
ing mass eigenstates. A detailed discussion is given in appendix [Al In general, the CP-odd 



4 



component S3 mixes with the CP-even fields S12 and the resulting mass eigenstates do not 
have a definite CP quantum number. If the scalar potential is CP symmetric this admixture 
disappears; in this particular case, A(x) = S3 (or) ancQ 



(3) 



Performing a phase redefinition of the neutral CP-even fields, we can fix the sign of sin a. In 
this work we adopt the conventions Mh < Mh and < a < it, so that sin a is positive. 



c: 


1= 


cos a 


sin a 






— sin a 


cos a 




2.1 Yukawa Alignment 

The most generic Yukawa Lagrangian with the SM fermionic content gives rise to FCNCs be- 
cause the fermionic couplings of the two scalar doublets cannot be simultaneously diagonalized 
in flavour space. The non-diagonal neutral couplings can be eliminated by requiring the align- 



ment in flavour space of the Yukawa matrices 25 ; i.e., the two Yukawa matrices coupling to 



a given type of right-handed fermions are assumed to be proportional to each other and can, 
therefore, be diagonalized simultaneously. The three proportionality parameters <j/ (/ = u, d, I) 
are arbitrary complex numbers and introduce new sources of CP violation. 

In terms of the fermion mass-eigenstate fields, the Yukawa interactions of the A2HDM 



read 25 



Y 



V2 

v 



H+ {u fa VM d V R - ^ MlVV L ] d + Q vM{P R l) 
-7,T,vfd [f M fV*f] + h - c -' 



(4) 



where Vr,l = ^r 5 - are the right-handed and left-handed chirality projectors, Mf the diagonal 
fermion mass matrices and the couplings of the neutral scalar fields are given by: 



VdJ 



Vu 



(5) 



1 In the usually adopted notation a — a — (3, where a is the rotation angle expressing the two mass eigenstates 
h and H in terms of the CP-even neutral fields of the original scalar basis 4>\{x) and ^(x). Since the choice of 
initial basis is arbitrary, the parameters a and (3 are in general unphysical; their values can be changed at will 
through SU(2) rotations. These angles only become meaningful in particular models where a specific basis is 
singled out (through a symmetry for instance). 
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Table 1: CP-conserving 2HDMs based on discrete Z2 symmetries. 



Model 






Q 


Type I 


cot [3 


cot /3 


cot /3 


Type II 


— tan f3 


cot 


— tan (3 


Type X 


cot (3 


cot /3 


— tan /3 


Type Y 


— tan j3 


cot /3 


cot (3 


Inert 












As in the SM, all scalar-fermion couplings are proportional to the corresponding fermion masses. 
This linear dependence on the fermion mass is characteristic of the A2HDM framework and does 
not hold in non-aligned 2HDMs with FCNCs. The only source of flavour-changing interactions 
is the Cabibbo-Kobayashi-Maskawa (CKM) quark mixing matrix V p4j. All possible freedom 
allowed by the alignment conditions is determined by the three family-universal complex pa- 



rameters <j/, which provide new sources of CP violation without tree-level FCNCs 25 . The 
usual models with natural flavour conservation, based on discrete Z 2 symmetries, are recovered 
for particular (real) values of the couplings as indicated in Table [I] 

Quantum corrections induce a misalignment of the Yukawa matrices, generating small FCNC 



effects suppressed by the corresponding loop factors 24 26 , 55 56 . However, the flavour sym- 
metries of the A2HDM tightly constraint the possible FCNC structures, keeping their effects 



well below the present experimental bounds 24 29 



The orthogonality of the rotation matrix TZ, implies the following relations among the 
Yukawa couplings of the three neutral scalars: 

33 3 



E (yf) 2 = 1 • E \yf\ 2 = 1 + 2 m 2 . E *f *« = 1 

i=l i=l i=l 
3 3 

Y y% n a = , Yl y ^ ni2 = ^ ' 



i=l i=l 
3 3 



Y y Z n ® = i w ' y ™ Uiz = ~ i ^- (6) 



i=l i=l 
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2.2 Bosonic Couplings 

The full set of interactions among the gauge and scalar bosons is given in appendix |Bj The 
relevant vertices for our analysis are the ones coupling a single neutral scalar with a pair of 



gauge bosons. As shown in Eq. (53), they are identical to their SM counterpart, with the field 
Si taking the role of the SM Higgs. Therefore (VV = W + W~, ZZ), 

9 v 9vv — T^n 9kvv > (7) 

which implies 

9hvv + 9hvv + 9 aw = (9hvv) ■ (8) 
The strength of the SM Higgs interaction is shared by the three 2HDM neutral bosons. In the 
CP-conserving limit, the CP-odd field decouples while the strength of the h and H interactions 
is governed by the corresponding cos a and sin a factors. Thus, a general feature of 2HDMs is 
that, at tree level, the couplings of the neutral scalars to vector bosons cannot be enhanced over 
the SM value and obey the custodial symmetry relation g^ zz — 9<p°ww- Observing a scalar 
boson with a somewhat enhanced coupling to vector bosons or a deviation from custodial 
symmetry [57] would therefore be in clear contradiction with the predictions of this class of 
models. The relations ^ and (pi) establish a connection between the couplings of the observed 
126 GeV resonance and searches for other neutral and charged scalars within the A2HDM. 

In order to compute the two-photon decay widths of the neutral scalars, one also needs 
their couplings to a pair of charged scalars, generated through the scalar potential discussed 
in appendix |A| Since these couplings depend on still unknown parameters, we will parametrize 
the corresponding interaction as 

CqPH+s- = -v 5^ \1h+h- ¥°i H + H~ . (9) 

Explicit expressions for the cubic couplings A^o H + H - , in terms of the Higgs potential parameters, 
can be found in appendix |A| If CP is assumed to be an exact symmetry, Xah+h- = 0. 



3 Higgs Signal Strengths 

The experimental data on Higgs searches is given in terms of the so-called signal strengths, 
measuring the observable cross sections in units of the corresponding SM expectations. At 
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the LHC, the relevant production mechanisms for a SM-like Higgs particle are gluon fusion 
(gg — > H), vector boson fusion (qq' — > qq'VV — > qq'H), associated production with a vector 
boson (qq' — > WH/ZH) and the associated production with a tt pair (qq/gg —> ttH). The 
Higgs decay channels explored so far are 77, ZZ^*\ WW^*\ bb and t + t~. 
In order to fit the experimental measurements, we consider the ratios : 

= a (PP <Pi) Br (^° 77) = -»■ JjV°) Br(^° -)> 77) 

^ 77 cr(pp -)> /i)sm Br(/i -> 77) S m ' <?{pp ->■ jjh) S M Br(/i ->■ 77)sm ' 

^0 _ o-(pp ->. ^) Br(y>9 -> VV) _ cr(pp -> jjy?°) Br(y?° ^ WW) 

^ = <r(PP->>OsMBr(fc-> WW ' = *(pp ^ jjfc) SM Br(fc^ WW) SM ' 

*? = °(PP ~+ Br(^° -» rr) ^ = -> V>°) Br(y° -» 66) 

^ <t(rp ^/i)sMBr(/i^rr)sM' ^ a(pp ^ Vh) S uBr(h ^ bb), ' 



where V = W, Z and j stands for jet. QCD corrections cancel to a large extend in these ratios, 
provided that a single production mechanism dominates. This certainly applies to p 7 j, jjy V and 
\i T T which are governed by the dominant production channel through gluon fusion. The same 
would be true for Pwwjj anc ^ Vyyjj (gauge-boson fusion), and /^y (associated production), 
assuming that there is no contamination from other channels. It is convenient to express the 
ratio of the branching fractions as: 

Br(y?° -» X) 1 r(y? -» X) 

Br(/i X) SM p(<P°) r(/i -> X) SM ' 1 1 

where p(<p°) measures the total decay width of the scalar 99° in units of the SM Higgs width, 

r(^°) = r S u(h) . (12) 

Particularizing to the A2HDM and assuming only one dominant production channel in each 
case Q one finds: 



i4v = (^i) 2 Mvd'T + MvdW p{^)~\ ^w„ = (^i) 4 p(^ 



0\-l ..<Pi . (<r> \4 „/,„0\-l 



fi\2 1 T™/„.^?\2a-2 



o^-l ..Vi . /^Vi \2 „/,„0\-l 



= c# c# p(^)- 1 , /i!;,, = (7^) 2 c% P (^r\ 



(13) 

The contamination of the different Higgs production mechanisms in h — > 77 (jj) is discussed in appendix C 



where (3f 



4m^/M 2 ) 1 / 2 . The one- loop functions are given by 



^99 





£ g Re(^)^) 


2 
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E 9 Im(^)/C(x 9 ) 2 


-> ^)SM 




2 



(14) 



and 



r(^° ->• 77) 

T(h -> 77)sm 



(15) 





2 

+ 








2 



with iV^ and the number of colours and the electric charge of the fermion f,Xf — Am 2 JM 2 



and xw = 4M^/M 2 ? . Notice that the ratios (10) are defined for M v q 



Mi 



/'SP, 



The two 



separate terms in the numerators of Eqs. (14) and (15) correspond to the CP-even and CP- 



odd structures ip^X ftv X ,lv and tp^X ttv X ,J ' v , with X^ v = G^ u (F^) in the gluon (photon) case 
and X^ v = e fll " 7p X ap . The functions J-(xf), JC(xf) and Q(xw) contain the triangular 1-loop 
contributions from fermions and W bosons. We will neglect the masses of the first two fermion 
generations. Since F(xf) and JC(xf) vanish for massless fermions, we only need to consider the 
top, bottom and tau contributions; the last two are negligible in the SM, but in the A2HDM 
could be enhanced by the alignment factors and q. In we have also considered the 
contribution from a charged-scalar loop parametrized by 

„,2 



2M 2 H± 



\°H+H~ A{XH1 



(16) 



with xh± = 4M^±/M 2 . The explicit expressions of the different loop functions are: 



JF(x) 
A{x) 



X 



[4 + (x-l)f(x)\ 



-x 



G{x) 
K{x) 



— 3x 
fix), 



-x 



x 2 ^jf(x) , 



(17) 



with 



/(*) 



-4arcsin 2 (l/^/a;) 



In 



— ITT 



X ^ 1 
X < 1 



(18) 
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4 Phenomenological Analysis 



We are interested in analyzing the current LHC and Tevatron data within the A2HDM. The 
experimental information on the new neutral boson is certainly in early stages; some decay 
channels have very big uncertainties while some others have not even been seen yet. Neverthe- 
less, while more precise information on all possible production and decay channels is necessary 
in order to make a detailed study, present data already allow us to extract significant constraints 
on the parameter space of the model. 

The deviations from the SM expectations originate from several sources. The three neutral 
scalars of the A2HDM have couplings to the gauge bosons which are different (smaller in 
absolute value) than the ones of the SM Higgs: in SM units they are given by IZn. The 
Yukawa couplings get also multiplied by the factors y^' , which are functions of TZij and the 
parameters qf. Moreover, the presence of a charged scalar manifests in one additional one- loop 
contribution to the (p® — > 27 decay amplitudes, parametrized through the constants C^ ± . In the 
limit of CP conservation, there are two clear candidates for the new scalar, the CP-even fields 
h and H (we will nevertheless analyze later the unlikely A possibility). The A2HDM allows 
in addition for physical CP-violating phases, both in the scalar potential and the Yukawa 
couplings, generating mixings among the three neutral scalars and CP-odd contributions to the 
Higgs-like signal strength parameters. Being quadratic in the CP-violating parameters, this 
last type of corrections could be expected to be small. However, the current bounds on the 



A2HDM couplings still allow for sizeable effects 24 29 



Sensitivity to the top-quark Yukawa coupling and to a lesser extent to the bottom coupling 
appears through the one-loop production mechanism of gluon fusion and in the 77 decay 
channel. Neutral scalar production via pp — > tLp®j(b) could provide complementary information 



on the top Yukawa coupling when more data becomes available 58 59 . The most important 
constraints on the bottom Yukawa coupling come indirectly from the total decay width, which is 
in general dominated by 9?° — > bb, and the measurement of scalar production with an associated 
vector boson (qq' — > (p®V — > (bb)V). Neutral boson production via top-quark fusion with 
subsequent decay into a pair of b quarks, qq/gg — > tiip® — > ti(bb), in which the bottom and top 
Yukawa couplings appear at tree level will also play an important role; the current experimental 



sensitivities in this channel are still low 60, 61 . The r Yukawa coupling is directly tested 
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through (p® — > t + t~ , the most accessible production mechanisms at the LHC being in this case 
vector-boson fusion, associated production with a vector boson and gluon fusion. 

For a given choice of neutral scalar-field candidate tp® and its couplings, we define the \ 2 
function as 

= E V ff 2 ^ > (19) 

where k runs over the different production/decay channels considered, /t^ and are the mea- 
sured Higgs signal strengths and their one-sigma errors, respectively, and ^ the corresponding 



theoretical predictions in terms of the A2HDM parameters, as given in Eqs. ( 10 ) and ( 13 ). Scan- 
ning over the allowed parameter space, we then look for those sets of couplings minimizing the 
X 2 and their corresponding uncertainties. The details about the statistical treatment and data 
used in this work are presented in appendix O 



We will first analyze the CP-conserving limit in section 4.1, where we will also study some 



particular scenarios often adopted in previous works. In section 4.2 we will discuss the most 
general case, without making any assumption about the scalar potential, and analyze the 
present constraints on the complex Yukawa couplings of the assumed 126 GeV scalar boson. 



4.1 The A2HDM in the CP-conserving limit 

Assuming that the Lagrangian preserves the CP symmetry, the two CP-even neutral scalars 
h and H couple to the gauge bosons with reduced couplings TZn = cos a and 72.21 = — sin a, 
respectively, and their Yukawa couplings are real: 

y h f = cos a + $f sin a , yf = — sin a + <r/ cos a . (20) 

The CP-odd boson A does not couple at tree-level to W + W~ and ZZ (72 3 i = 0), while its 
fermionic couplings are purely imaginary (pseudoscalar interaction): 

Vd,i = ifyi, Vu = (21) 

4.1.1 A light CP-even Higgs at 126 GeV 

We will first focus in the most plausible possibility that the lightest scalar h corresponds to 
the observed neutral boson with = 126 GeV. The alternative choice of the heavier field H 
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can be easily recovered through an appropriate change of the mixing angle, a — > a — tt/2, and 



will be further discussed in section 4.1.5 We will also consider later, in section 4.1.6, the more 



exotic case of a CP-odd Higgs A. In this first analysis we assume that the charged scalar is 
either very heavy or its coupling to the neutral Higgs is very small, so that its contribution C^ ± 
to the h — y 77 decay width is negligible. We also assume that the bounds from flavour physics 
are naturally evaded, as it is the case at large values of the charged scalar mass. The H ± 
contribution to the diphoton decay width as well as the flavour constraints will be considered 
later in section 14.1.31 

The minimization of x 2 (/i) leads to two different solutions, differing in the sign of the 
top Yukawa coupling. The central values of the corresponding A2HDM parameters and their 
statistical one-sigma errors obtained from the global fit are: 

cosa = 0.99t°;°6> y h u = 0.8i°i , = 0.7± 0.3 , |yf | = 0.8 ± 0.5 , (22) 

and 

cos 5 = 0.99^, y5 = -0.8igi, \y h d \ = 1.1 ±0.3, \y* | = 0.9 ± 0.5 . (23) 

In both cases, the gauge coupling ghvv is ver Y close to the SM one. Changing simultaneously 
the signs of cos a and yj leads obviously to identical Higgs signal strengths and, therefore, to 
two equivalent solutions. 

In the first solution the W and top-quark loops contribute with different signs to the 
h —> 77 amplitude, giving a destructive interference as in the SM. The needed enhancement 
of the 27 branching ratio is obtained through a smaller total decay width, p(h) w 0.6. This 
pushes upward the ratios /A, and /jtt^j, allowing to explain part of the excess experimentally 
observed in these two channels. However, the gluon-fusion production channel has a smaller 
cross section than in the SM. The combined effect results in a small increase of the 77 channel, 
ps 1.1, while a much larger enhancement remains in the 77 jj case, fJ>!t~jj ~ 1.5. 

The second solution corresponds to a top-quark contribution to h — > 77 with the opposite 
sign, so that it interferes constructively with the W amplitude. This allows one to explain the 
27 excess without hardly modifying the total decay rate, p(h) ps 1.1 and providing a slightly 
better fit. 
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Figure 1: Global fit to the A2HDM, in the CP-conserving case, in the planes y^ — y% (left) and 
Vu ~ 2/f (right). The parameters not shown in each case are fixed to the best global-fit point. The 
orange, yellow and gray areas denote 68%, 90% and 99% CL regions. The dashed lines correspond to 
fixed values of ffi . 

In both solutions there is a sign degeneracy in the bottom and tau Yukawa couplings. 
Although the tree-level decays h —¥bb and h — > t + t~ are insensitive to these signs, the loop- 
induced processes gg — )■ h and h — > 77 receive contributions from the bottom and tau (only the 
77 process) Yukawas, which interfere with the leading top and W ± (in the 77 decay) amplitudes 



as shown in (14) and (15). In the SM the bottom and tau contributions are negligible, but their 



effect could be relevant in the A2HDM if the top Yukawa coupling is considerably suppressed or 



if the parameters $d,i are large. However, this is not the case for the fitted Yukawa values in (22) 



and (23), which are of 0(1) for both solutions, leaving the sign of the bottom and tau Yukawas 
undetermined. The relevance of the t + t~ and bb channels to determine possible deviations 
from the SM and within the different Z 2 versions of the 2HDM, which could be pointing to a 
more general Yukawa structure as provided by the A2HDM, has been emphasized recently in 
Ref. [52]. 

In Fig. [T]we show graphically the results of this global fit, giving the allowed regions in the 
Vu ~ Vd (left) an d Vu — Di (right) planes at 68%, 90% and 99% CL. The parameters that are 
not shown are, in each case, set to the best global-fit point. The sign degeneracy in the r and 
b Yukawa couplings is clearly observed. Moreover, the right panel shows a somewhat reduced 
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Figure 2: Allowed ranges for the Higgs signal strengths obtained from the fit (23) at la (black, dark) 
and 2a (blue, dark), together with the averaged experimental data from the ATLAS, CMS, CDF and 
D0 collaborations with the corresponding la errors (orange, light). 

sensitivity to the leptonic coupling y\ . The SM-like solution {y^Vdi) = (1) 1) ues inside the 
90% CL allowed region; however, at 68% CL the top Yukawa has the sign flipped with respect 



to the SM, i.e., only the solution (23) remains. Similar results have also been obtained in 



Ref. 52 



The allowed ranges, at the la and 2a level, for the different Higgs signal strengths in the 



fit (23) are compared in Fig. 2 with the experimental values. A good agreement with data is 



obtained in all cases. Previous analyses within the CP-conserving A2HDM have been performed 



in Refs. 52 , 53 , using a different notation, also finding good agreement with the data. 



Using the sum rules in Eqs. ^ and (|8]), we can extract constraints on the heavy CP-even 
Higgs couplings from our global fit with = 126 GeV. For the solution (23 ) we find at 68% CL 



that the coupling of H to vector bosons is suppressed, sin a < 0.37, while its coupling to top 
quarks is very large, \y^\ > 4.6. This region of parameter space requires a very large value of 
\q u \ in order to flip the sign of y\, which is the top Yukawa of h. Such large values of |sj would 
then imply a significant enhancement of the production of H via gluon fusion and can give rise 
to non-perturbative H + ib, Hit and Ait couplings. This was noted previously within the same 



context in Ref. 52 . 
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Figure 3: Global fit within 2HDMs of types I (upper left), II (upper right), X (lower left) and Y 
(lower right), at 68% (orange), 90% (yellow) and 99% (gray) CL. The dashed lines correspond to 
constant values of /x^ 7 . 

4.1.2 Global fit within Z 2 models 

The usual 2HDMs with natural flavour conservation, based on discrete Z 2 symmetries, are 
particular cases of the CP-conserving A2HDM, with taking the values given in Table [Tj 
Thus, the three alignment factors are determined by a single parameter through the constraints 
W = ?d = ?j = cot/3 (type I), ^ = = -cf 1 = cot/3 (type II), q u = q d = -q- 1 = cot/3 



(type X) and 



-i 



q = cot/3 (type Y), with cot/3 = Vi/v 2 > 0. This leads to specific 



relations among the production cross sections and decay rates for the Higgs bosons that can 
be tested with the LHC data. The separate measurement of the various Higgs signal strengths 
should allow to disentangle the different scalings of the three Yukawa couplings. In particular, 
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Figure 4: Allowed ranges for the Higgs signal strengths in 2HDMs of type I, II, X and Y, at la 
(black, dark) and 2a (blue, dark). Other captions as in Fig. [1| 

exclusive Higgs production measurements in the final states r + r~ and bb will be crucial to test 



the different Z 2 versions of the 2HDM 32,42,52 



Figure [3] shows the results of the global fit for the 2HDMs of types I, II, X and Y, assuming 
that the lightest neutral Higgs h is the boson observed around 126 GeV. Allowed regions at 
68%, 90% and 99% CL are shown, together with lines of constant ^ . The relevance of the 
diphoton channel is evident from the figure. In models I and X, an allowed region around 
cos« ~ 1 appears, where there is no sensitivity to q u since its contribution to the neutral 
Yukawa couplings is suppressed by sin a; in this region the couplings of h to vector bosons and 
fermions are close to the SM ones. Another allowed region appears for negative values of cos a, 
in which the W ± and top-quark loops contribute with the same sign to the h — > 77 decay 
amplitude, thus allowing for a constructive interference. Both solutions with cos a ~ ±1 are 
present for the inert model (type I with <^ u = 0). There is a third allowed region at large values 
of the top Yukawa and negative cos a, which approaches cos a = — 1 as <; u increases. 

In models II and Y the solutions around cosd ~ ±1 reduce to two extremely narrow vertical 
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lines and one small region at low q u and positive cos a, which remain allowed at 99% CL but 
are not present at 90%. The solution at large values of the top Yukawa and negative cos a is 
also present, but in a region much smaller than in models I and X. 

Figure [4] shows the allowed ranges for the Higgs signal strengths obtained in these four types 
of 2HDMs (I, II, X and Y). The agreement with the data is good; however, as already noted in 



Ref. 42 , the preferred region has large values of \q u \, which are ruled out from flavour physics 
constraints for a charged Higgs boson below the TeV scale. Large values of \<; u \ can also make 
some top-quark Yukawa couplings non-perturbative, as commented in the previous section. 

4.1.3 A charged Higgs and the diphoton excess 

One of the most distinctive features of 2HDMs with respect to other alternative scenarios of 
electroweak symmetry breaking is the presence of a charged scalar boson in the spectrum. The 
present experimental lower bound on the H ± mass is Mjj± > 80 GeV [62] (assuming that the 
charged scalar H + only decays into fermions). Direct searches for charged Higgs bosons at the 



Tevatron 63 and the LHC [64] have also been performed with null results so far. Current LHC 
data are sensitive to such charged scalar through the h — > 77 decay channel. The one-loop H ± 
contribution can interfere with the W and fermionic amplitudes, thus being able to enhance 
or suppress the decay rate. The exact value of the charged Higgs contribution C^± depends 
on the cubic Higgs coupling \hH+H- and the charged Higgs mass M H ±. One expects however 
that \Cx± \ < 0(1) based on perturbativity arguments (see appendix [d|). 

When considering a relatively light charged Higgs boson, one must take into account con- 
straints from electroweak precision tests and the flavour sector; a light H would contribute 
sizably to loop-induced processes, such as Z — > bb, b — > 57 or B°-B° mixing. These phe- 
nomenological constraints have been analyzed in detail within the framework of the A2HDM 



in Refs. 26-29 , where it has been found that a charged Higgs below the TeV scale would 



require \s u \ < 2 to be compatible with present data. This rules out the hypothetical scenario of 



a top Yukawa coupling with flipped sign, as found in (23) and also favoured by the fits shown 
in Fig. [3] within the four types of Z 2 models. The reason is that current h — > WW, ZZ, "fy(jj) 
data require | cos<5| ~ 1 (i.e., the gauge coupling of the new neutral scalar should be close to 



the SM one). Since the top Yukawa coupling is given by = cos a + <, u sin a, in order to flip 
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the sign of y% one needs then a large value for \q u \, which is excluded by the previous bound. 

Including the charged-Higgs contribution, it is no longer necessary to flip the sign of the 
top Yukawa in order to enhance the h — > 77 decay width. The best fit region is now obtained 
for Yukawa and gauge couplings close to the SM limit: 



cos a 



0.98 



+0.02 

-0.06 > 



C h H± = (-2.8 ±1.3) U (16.0 ±1.3) 



1.0 ±0.2. 



\y*\ 



1.1 ±0.3 



\Vi 



0.8 ±0.5. 



(24) 



The two disjoint C^ ± solutions correspond to either a constructive interference of the H and 
W ± amplitudes or a destructive one but with a charged-Higgs contribution so large that it 
reverses the sign of the total h — > 27 amplitude. In both cases, one obtains a better fit than 



in the SM and also better than the previous A2HDM fits (except for (23) which is comparable 
to this one). The presence of the charged Higgs allows one to easily explain the h — > 77 (jj) 
excess without large modifications of the total decay rate (i.e., p(h) « 1.1). The fit predictions 
for the fik ratios and their one and two-sigma statistical errors are shown in Fig. |5j In all cases, 
good agreement with the data is obtained. 

In Fig. [6] we show the allowed regions of the (\XhH+H- \ ,M H ±) plane, corresponding to the 
two possible fitted values of C^ ± , at 68% and 90% CL, together with the perturbativity bounds 
discussed in appendix [Dj Clearly, the solution with a very large contribution to h — > 77 from the 
charged Higgs (C^± ~ 16) is excluded if one requires the theory to be perturbative. We obtain 
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Figure 6: Allowed regions of the (\^hH+H- \ )-^ff±) plane, corresponding to the two possible fitted 
values of Cjt±, at 68% (orange, dark) and 90% CL (yellow, light). The blue (hashed) area, between 
the left vertical axis and the dashed line, is the domain where the theory remains perturbative. 

an upper bound for the mass of the charged Higgs around 300 GeV, at the one-sigma level. 
However, the bound disappears at the two-sigma level because the charged-Higgs contribution 
becomes compatible with zero. 



Perturbativity 




4.1.4 Inert 2HDM 



In the inert 2HDM a Z 2 symmetry is imposed, in the Higgs basis ([2|, under which all SM 
fields and $i are even while $2 — > —&2- Terms with an odd number of $2 fields in the scalar 



potential (30) are then forbidden by the Z% symmetry, therefore ^ = Xq = A 7 = 0. In this 
case there is no mixing between the CP-even neutral states h and H, and the scalars H, A 
and H ± decouple from the fermions. The couplings of the remaining Higgs field h to fermions 
and to vector bosons are the same than in the SM (i.e., cos a = 1 and = 1). Thus, only 
the diphoton channels can show a deviation from the SM prediction (assuming that there are 
no open decay channels other than the SM ones). From the global fit of this scenario, we 
find a charged-Higgs contribution to the h — > 77 amplitude in the range C^± G [—1.7, —0.89] 
at 68% CL and C h H± G [-2.4,-0.1] at 90% CL. We have assumed that M H ± is greater than 
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Mh/2 ~ 63 GeV so that C^ ± is real; for lower charged- Higgs masses, it would develop and 
imaginary absorptive part. The fitted negative sign of C^± causes a constructive interference 
with the W ± amplitude in the h — > 77 decay width. 

Note that in the limit = 0, the charged Higgs does not couple to fermions independently of 
any assumption on the scalar potential, see Eq. Q. The implications of this more general case 
for the neutral Higgs boson phenomenology as well as the possibility of a very light charged 
Higgs boson are discussed in section 4.2.2| Detailed analyses of the inert 2HDM and the 



possibility of a Dark Matter candidate within this model, in light of the LHC data, can be 



found in Refs. 38,65 



4.1.5 A heavy CP-even Higgs at 126 GeV 

We have discussed so far the phenomenology of the lightest Higgs boson, but there is nothing a 
priori preventing the boson discovered by ATLAS and CMS to be identified with the heaviest 
CP-even state H or with the CP-odd Higgs A. These possibilities have been already discussed 



in Refs. 34 36 52 . An analysis in terms of the more general CP- violating scalar potential, 
setting limits on the scalar-pseudoscalar mixing, has been done in Ref. |35| . 

Using the previous fits for h, it is straightforward to analyze the possibility of having a heavy 
Higgs with Mh = 126 GeV. Assuming that non-SM decays like H — >■ hh are kinematically 
forbidden or very suppressed, the constraints on the heavy Higgs boson couplings can be easily 
obtained from those of h through an appropriate change of the mixing angle: a — > a — tt/2. In 
this case the coupling of the heavy Higgs to vector bosons is close to the SM limit (sin a ~ 1), 
while the light-scalar ghvv couplings are suppressed by cos a ~ 0. The absolute values of the 
Yukawa couplings and all the other parameters remain unchanged. 

The LEP searches for neutral Higgs particles could have missed the light scalar h, since the 
associated production with a vector boson would be strongly suppressed. Moreover, \y^\ ~ |^| 
could be small enough to avoid the constraints from the usual h — > bb search mode. The OPAL 
collaboration performed a decay-mode-independent search for a light neutral scalar and found 
upper limits for the Higgs-strahlung cross section in units of the SM: (Jin) 2 = {dhvv / 9hvv) 2 < 



0.1 for M h < 19 GeV, and (^n) 2 < 1 for M h < 81 GeV 66 . Together with the constraints from 



electroweak precision tests at the Z peak, this provides useful information on the allowed mass 
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Figure 7: Left-panel: Constraint in the (M H ±, Ma) plane from the oblique parameters S, T and U. 
Right-panel: Constraints from the invisible Higgs decay width in the (\\Hhh\ , Mh) plane, assuming 
SM couplings of H to fermions and vector bosons. The orange (dark) and yellow (light) regions are 
allowed at 68% and 90% CL. 

spectrum for the remaining scalars. Using the current bounds from the oblique parameters S, 



T and U 67,68] (the corresponding A2HDM formulae are given in appendix |Ej), we show in the 

left panel of Fig. [7] the allowed regions in the (M H ±, Ma) plane. We have set M# = 126 GeV 

and sin a G [0.7, 1]. The constraints shown in the figure turn out to be determined by the T 

parameter, since S and U give weaker restrictions. 

A light neutral boson h or A below the kinematical threshold of Mff/2 ~ 63 GeV would 

have important phenomenological consequences, because the 126 GeV Higgs could decay into 

lighter scalars. These decay channels can be included in our fit in terms of an invisible decay 

width as long as we neglect possible contributions from cascade decays into the observed final 

states]^] In general one would expect in this case a suppression of the measured Higgs decay 

rates compared with the SM, due to the larger total width of the scalar H. Current data for 

the 77 channel, however, shows a slight enhancement over the SM prediction, thus placing 

strong bounds on possible invisible decays of the 126 GeV Higgs boson. Assuming that the 

heavy-Higgs couplings to fermions and vector bosons are SM-like (i.e., = 1 and sin a = 1), 

the best fit point is obtained for a null invisible H decay width; at 68% CL (90% CL) we obtain 

3 These effects are beyond the scope of the present work, but they could be relevant. For example, H —> 
A A — > 77 + 77 could be mistaken by a two-photon signal when the photon pairs are very collimated |69 
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an upper bound of 9% (20%) on the invisible H decay width (in units of the SM total decay 
width). 

Considering the scenario of a very light CP-even Higgs h, the decay width of the heavier 
CP-even scalar into hh is given by 

v 2 \* / 4M^ 1/2 



T(H hh) = — 1 f , (25) 

where the cubic scalar coupling Xuhh is expressed in units of v and can be obtained from 
Eq. (42). In the right panel of Fig. [7] we show the constraints from our T(H — > hh) fit in the 



(\^Hhh\, Mh) plane. Strong bounds are obtained for the cubic Higgs coupling, |Ajjh/i| ^ 10 -2 , 
as expected. 

Recent updates from the ATLAS collaboration in the high-resolution channels report a 
significant difference in the mass of the neutral boson as determined from H — > ZZ^ — >■ A£ 
(123.5 ±0.8 ±0.3 GeV) and H -»■ 77 (126.6 ±0.3 ±0.7 GeV) events (l). Here we do not consider 
as a possible explanation for this discrepancy, the possibility of having two quasi-degenerate 
Higgs bosons, since the current mass value in the H — » ZZ^ — > 4£ channel obtained by CMS, 
126.2 ± 0.6 ± 0.2 GeV [3], does not support this hypothesis. 

4.1.6 Degenerate CP-even and CP-odd Higgs bosons at 126 GeV 

A CP-odd scalar does not couple at tree level to two vector bosons; its decay to gauge bosons 
starts at the one-loop level and it is therefore very suppressed. For this reason, a pure CP-odd 
Higgs boson is already strongly disfavoured by present data as a candidate for the 126 GeV 
boson. However, the observed signal could result from two Higgs bosons with quasi-degenerate 
masses; this could explain the excess of 77 events observed by ATLAS and CMS. This possibility 
was proposed in Ref. [70] within the non-minimal supersymmetric extension of the SM, and 



has also been considered within the context of 2HDMs, both for Z2 versions 71 - 73 and with 



a more general Yukawa structure [52] . Model- independent methods to test experimentally for 



such possibility have also been proposed recently in Refs. 74 , 75 . 

We consider in this section the possibility of two Higgs bosons with quasi-degenerate masses 
around 126 GeV, one of them being CP-even and the other one CP-odd. We perform a global 
fit of the data with = M4 f=s 126 GeV, and comment on the alternative possibility of quasi- 
degenerate H and A. The observed Higgs signals strengths will then receive contributions from 
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Figure 8: Constraints in the (M H ±, Mh) plane for the case Mh = Ma = 126 GeV (left) and in the 
(Mjj±, Mh) plane for the case Mh = Ma = 126 GeV (right), from the oblique parameters S, T and 
U. The orange (dark) and yellow (light) regions are allowed at 68% and 90% CL. 

both particles: 

^ +A) = + ■ (26) 

Given the presently large experimental uncertainties, we neglect the small AVV coupling gen- 
erated at one loop. Therefore, among all the channels considered in this work, the CP-odd 
Higgs will only contribute to A — )■ rr and A — > 77. In both cases the dominant production 
channel is the gluon-fusion one. The loop-induced decay A — > 77 is only mediated by fermions. 
In Fig. [8] (left) we show the constraints on M H ± and Mjj obtained from the oblique parame- 
ters. These masses are varied in the ranges Mn± 6 [10, 10 3 ] GeV and Mh 6 [126, 10 3 ] GeV, 
while the coupling of h to vector bosons is kept close to the SM limit (i.e., | cosa| G [0.8, 1]), 
as suggested by the current experimental data. In the right panel of Fig. [8] we show similar 
bounds on the plane (M H ±, Mh), keeping the light scalar mass below Mh = Ma = 126 GeV 
and taking sin a G [0.8, 1]; in this case the oblique parameters require the existence of a charged 
Higgs below the electroweak symmetry breaking scale v = 246 GeV. 

In the scenario M^ = Ma = 126 GeV, the best fit region in the A2HDM parameter space 
is given by: 

cos a = 0.98 ±0.2, q u = -lA^H, |%| = 1.2 ± 1.2 , Q = -0.2ig;|. (27) 

The corresponding allowed ranges for the Higgs signal strengths, at la and 2a, are shown in 
Fig. |9j We obtain a smaller total decay width of the CP-even boson, p(h) ~ 0.7, which produces 
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Figure 9: Allowed ranges for the Higgs signal strengths from the global fit within the CP-conserving 
A2HDM for the case of degenerate Higgs bosons with Mh = Ma = 126 GeV, at la (black, dark) and 
2<T (blue, dark). Other captions as in Figure^ 

a sizeable enhancement of the A^ 7 jj signal strength (the CP-odd boson A does not contribute to 
this channel). On the other hand, the excess in the two photon channel comes from the decays 
of both A and h, which give contributions of similar size (/A, ~ /j,^ ~ 0.7). The remaining 
contribution of A is to the r + r~ decay channel, which is small (q is small). We must also notice 
that solutions with a flipped relative sign between the W and top contributions to h — > 77 
are not allowed because they would require large values of q u ; this would increase C^ g and 
generating a large excess in the r + r _ and 77 channels, exceeding the current experimental 
bounds. 

4.2 The CP-violating A2HDM 

In the A2HDM the up and down-quark as well as the leptonic Yukawa couplings are all indepen- 
dent complex parameters. Thus, one can expect a very rich phenomenology associated to the 
Higgs sector responsible for the breaking of the electroweak symmetry. Moreover, if one consid- 
ers the most general scalar potential, the neutral scalars h, H and A are not CP eigenstates but 
rather a mixture of CP-even and CP-odd fields, parametrized by the general orthogonal matrix 
1Z introduced in section [2] Thus, there are new sources of CP violation, both from the Yukawa 
sector and the scalar potential, which could lead to interesting phenomenological predictions. 
The study of CP-violating observables is beyond the scope of the present work and we 
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will defer it to future publications^] Nevertheless, we shall investigate next, the sensitivity 
of the different (CP-conserving) Higgs signal strengths to the CP-violating phases. Since the 
present data are consistent with the SM within rather large uncertainties, we will consider 
separately the different CP-odd possibilities, by fitting some complex coupling constants to the 
Higgs-signal-strength data while setting the remaining parameters to their SM-like values. 

4.2.1 Complex Yukawa couplings 

Let us consider ip® to be the observed boson with a mass of 126 GeV. We will analyze three 
simple scenarios that will serve to determine the sensitivity to its complex Yukawa couplings and 

to what extent the SM limit is preferred by present data. We will set two Yukawa couplings to 

ip° 

their SM values (yj* = 1), and find the preferred values for the remaining Yukawa coupling by 



minimizing the y 2 function. Figure 10 shows the resulting allowed regions for the top, bottom 
and tau Yukawa couplings when the coupling of Lp\ to vector bosons is fixed to 7£ji = 0.95; this 
value lies well within the 90% CL allowed band obtained from our previous fits. 

Since all the observables considered are CP-even, the bounds obtained are symmetric under 
Im(yp) — > — Im(?/J*). Moreover, the real and imaginary parts of the Yukawa couplings do not 
interfere. The sensitivity to Im(?/^) is similar to that obtained previously, when considering 



only real couplings. For tree-level decays this is obvious from Eq. (13), given that the parameter 
f3f is very close to one for f = b,r. For loop-induced decays this can be understood by observing 
that the loop functions ( p~T| ) are closely related, T(t) = 2r + ^/(r) + /C(r). For b quarks and 
r leptons, J~{Tf) ~ /C(ry); for the top quark there is a small but sizable difference between 
the contributions of its real and imaginary Yukawa parts. Note that in the limit IZn = 1 the 
Yukawa couplings of (p® become SM-like (yj' = 1) due to the orthogonality of 7Z; thus, there is 
no sensitivity to the parameters when considering the neutral Higgs couplings. The charged 
Higgs couplings on the other hand are proportional to ^/ and do not depend on the mixing 
matrix TZ. 



In the left upper panel of Fig. 10 we show the results of the fit for a complex top Yukawa 



coupling, while setting y^ = yf* = 1. The dashed lines show contours of constant value 



4 For theoretical studies about the CP-properties of extended Higgs sectors at the LHC and in possible future 



colliders see Ref. 76 and references therein. 
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Figure 10: Allowed regions at 68% (orange), 90% (yellow) and 99% (grey) CL for the complex top 
(upper-left), bottom (upper-right) and tau (lower) Yukawa couplings. In each plot the two Yukawa 
couplings not shown are set to their SM value and the coupling to vector bosons is taken to be IZn = 

ufi if 

0.95. The dashed lines show contours of constant values for (top plot), fJ* b £y (bottom plot) and 

<p° , , , 
fi T l rV (tau plot). 

for fj^if. The SM-like point (Re^*), Im^*)) = (1,0) lies outside the 90% CL region, but 
becomes allowed at 99% CL. It can be seen that the allowed region at 90% CL accommodates 
an enhanced 77 rate between one and two times that of the SM. Within this 90% CL region, 
p(<Pi) = 1.00 ±0.03 as expected, since the dominant decay channel is bb; the gluon fusion cross 
section is slightly reduced compared with the SM [C g g = 0.87 ± 0.28), while the 77 partial 
decay width is enhanced (C 7 4 = 1.67 ± 0.56). The preferred allowed region is that for which 
the top Yukawa coupling has opposite sign to lZn, thus, creating a constructive interference 
with the vector boson contribution for the ip® — > 77 amplitude. The other option would be to 
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in 

have a significant imaginary component lm[y u i ), which would also enhance the 77 rate. 



The right upper panel of Fig. [10] shows the fitted values for the complex bottom coupling, 
with the top and tau Yukawa couplings set to their SM values. The dashed lines indicate 
contours of constant value for fi^y In this case the SM limit (Re^*), Im^y/)) = (1,0) lies 

in" 

inside the 90% CL allowed region, which accommodates 0.7 < n b l v < 1.2. In this 90% CL 
region, the total decay width is rescaled by p((p®) = 1.11 ± 0.67; the gluon-fusion cross section 
ratio is C g g = 1.15 ± 0.10, while the 77 partial decay width turns out to be slightly suppressed 
with respect to the SM, C# = 0.89 ± 0.10. Since the total decay width depends strongly on 
the value of | 2 , a large variation range is obtained for p(<£>°). 



In the lower panel of Fig. [LOj we show the fitted values of the complex r Yukawa coupling 
assuming y^ = y^ = 1. Contours of constant value for fi^y are also shown as dashed lines. 
We obtain that the signal strength fM^y < 1.5 lies within the 68% CL allowed region. The total 
Higgs decay width and the gluon-fusion cross section are equal in this case to the SM ones, 
while some suppression is observed in the 77 partial decay width: at 90% CL, C 7 4 = 0.90±0.11 
is obtained. This scenario is therefore disfavoured by the observed excess in the two-photon 
channel. 



4.2.2 A fermiophobic charged Higgs 

In the limit — > the charged Higgs does not couple to fermions, independently of any 
assumption about the scalar potential. Such fermiophobic charged Higgs could have avoided 
detection at LEP while being very light. Current LHC searches, as well as searches at the 
Tevatron, would have also missed such particle since it can neither be produced via top decay 
nor decay into fermions. Flavour constraints on this charged Higgs are also avoided trivially. 
Detecting such particle in an experiment is therefore quite challenging, since it can only be 
produced in processes involving vector bosons and/or neutral Higgs particles; the same occurs 
for its decay channels. 

The case of a fermiophobic charged Higgs is however highly predictive in the neutral Higgs 
sector, since all the channels which do not involve the 77 (jZ) final state only depend on 
one free parameter, Tin. The rescaling of the Higgs coupling to vector bosons in this case is 
the same as that of the neutral Yukawa ones, y^ = g^oyv 1 9 S ^vv = which implies that 
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all Higgs signal strengths are rescaled by a factor IZ^ with respect to the SM, meaning that 
/•*&& = V&r = Hwwzz = Pi.vi)~ l ^-ii = ^li> m an y °f the relevant production mechanisms. 
Therefore, in this scenario the signal strengths of the three neutral scalars are correlated: 

^// = ^ww,zz = !• (28) 

<p?=h,H,A ip°=h,H,A 

Present data on the neutral Higgs boson are sensitive to a fermiophobic charged Higgs 
through the loop-induced decay (p® — > 77. The charged-scalar contribution to this decay can 
be sizeable for a light H ± , and this is a quite interesting situation in view of the possibility 
to detect such particle in the future. Assuming that the scalar with a mass of 126 GeV does 



not decay into lighter scalars, we show in Fig. 11 the allowed region in the parameter space 
(7£;i,C^ ± ). For the x 2 fit we have only considered real values of C^±, which is true above 
the kinematical threshold M#± > M m o/2 63 GeV, as we have mentioned before. In the 
figure we also show dashed contour lines of constant /x 77 . It can be observed that the preferred 
relative sign between the charged Higgs and the W ± contributions to the 77 decay rate is such 
that it causes a constructive interference, thus enhancing slightly the 77 decay rate. The fit 
prefers a gauge coupling close to the SM one (Xmin * s obtained for Tin ~ 0.95) and puts the 
90% CL lower bound \K a \ > 0.79. The SM-like point (TZn,C^ ± ) = (1,0) lies outside the 68% 
CL region, but is allowed at 90% CL (although close to the boundary). The presence of a 
non-zero (and negative) Cjj± contribution is clearly favoured, while the preference for a slightly 
reduced gauge coupling implies a small suppression of the total decay width compared with the 
SM {i.e., p(ip?) = 0.85 ± 0.19, at 90% CL). From the global fit, $\ = = 1.45 ± 0.49 is 

obtained at 90% CL; all the other Higgs signal strengths that are not affected by the charged 
Higgs contribution are equal to fi = 0.8 ± 0.2. 



The right panel in Fig. 11 shows the corresponding allowed regions in terms of the variables 
\ w o H+H - and M H ±. The value of IZn has been set to its best fit point. Also shown in the 
figure, is the region satisfying the perturbativity constraints discussed in appendix |D} 

For the previous discussion we have not made any assumptions on the quantum numbers 
of the scalar field <£>°; we have only assumed that M^o = 126 GeV and that its decay into 
lighter scalars is not allowed. Thus, the obtained results are general and apply both to a CP- 
conserving and to a CP-violating scalar potential. It must be noted that in the limit \7Zn \ = 1 
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Fermiophobic charged Higgs Fermiophobic charged Higgs 




^il l^hH + i/-| 

Figure 11: Allowed regions at 68% (orange), 90% (yellow) and 99% CL (gray) for a fermiophobic 

in ip° 

charged Higgs on the parameter space (lZn,C^ ± ); dashed lines denote contours of constant /j,^ (left). 
The right plot shows the corresponding 68% and 90% CL regions in the parameters ^ v o H + H - and 
Mjj± , setting the value of TZn at its best fit point. The region where perturbation theory remains valid 
is indicated in blue (hashed). 

the phenomenology of (p® becomes identical to that of the SM in every channel, except for 77 
and 7Z which are affected by the if ± contribution. For a fermiophobic charged Higgs lighter 
than M v o/2 « 63 GeV, C£ ± develops an imaginary absorptive part. If kinematically open, 
the channel (p® — > H + H~ would increase the total width of the Higgs boson; furthermore, in 
this scenario the production cross section is always less or equal to the SM. Therefore, the 
signal strengths would be reduced in every channel, with respect to the SM. This is in clear 
contradiction with the data, specially with the measurements for the two-photon channel. 



4.2.3 CP-even and CP-odd neutral scalar mixing 

A CP-violating scalar potential generates mixings among the three neutral scalars, which are 
no longer CP eigenstates. Here, we are interested in exploring the possibility that the observed 
126 GeV state could be the CP-odd scalar with a small CP admixture of the CP-even ones. A 
similar analysis within 2HDMs of types I and II, with explicit CP violation and soft breaking of 



the Z2 symmetry has been done in Ref. 35 , placing numerical bounds on the size of a possible 



CP-odd component for the scalar particle with 126 GeV of mass. 
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In the presence of CP violation, the admixture between the three neutral scalar fields is 
described by the 3-dimensional orthogonal matrix 1Z which diagonalizes their mass matrix. 
This diagonalization can be done numerically, once the parameters of the scalar potential are 
known, but a simple analytical solution is not available for the most general case. It is well 
known, on the other hand, that in the CP-conserving limit the mass-matrix simplifies and it 
is possible to give explicit expressions for the masses and physical states in terms of the scalar 
potential parameters. A reasonable assumption when dealing with the general 2HDM scalar 
potential, is that the CP-violating terms are small; this makes a perturbative expansion in these 
parameters a valid approximation in principle. In appendix [A] we provide explicit analytical 
expressions for the neutral scalar masses and the corresponding eigenstates to leading order 
in the CP- violating parameters of the scalar potential The corrections to the masses are 
quadratic in Agg, while the mixing between the CP-even and CP-odd states is only suppressed 
by one power of X\ e , making this effect the dominant one. 

Let us assume that the discovered boson is the state A = S3 + H31S1 + H32S2, with TZ^i 
and 7^.32 the small CP-even admixture coefficients. To simplify the discussion, we consider a 
simple scenario in which we set the parameters ^ U)( y = 0. The Yukawa couplings, as well as the 
coupling to vector bosons, are equal in this case, = IZ31. From a global fit to the data, we 
find a lower bound on the admixture coefficient: 72.31 > 0.83, at 99% CL. This result is mainly 
driven by the measurements in the VT + W /_ , ZZ and 77 channels, which are SM-like to a good 
degree. 

We can analyze whether such large values for the correction 72-31 can be obtained for natural 



values of the scalar potential parameters. From Eq. (47), one has: 



v 



31 ~ {M\ - Ml) {Ml -Ml)' [ J) 

Thus, large mass differences between the scalar states suppress the effect of mixing due to 
CP violation in the scalar potential; on the other hand if the scalar bosons have very sim- 
ilar masses these effects could be considerably enhanced. Assuming that {X^qI < 10 _1 we 
obtain 72 3 i < [(Ml -Ml) (Ml - Ml)]' 1 10 8 GeV 4 , which implies that |72. 3 i| < lO" 2 for 
M H > M h > 300 GeV. Of course, when either M h ~ M A or M H ~ M A the coefficient 
72.31 diverges and the approximations used in appendix [A] are no longer valid. The general 
formalism to describe the dynamics of CP violation near degenerate neutral Higgs bosons has 
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Figure 12: Allowed regions at 90% CL (yellow) on the parameter space (Hn, "^3), for real alignment 
parameters in the ranges \q u \ < 2 and \$d,i\ < 10 (left). The right plot shows the corresponding 90% CL 
region for the parameters (Tli2, Ufa)- 



been developed in Refs. 77, 78 . In Ref. [79] the effect of resonant enhancement of H and A 
mixing was studied for the CP- violating 2HDM in the decoupling limit, M\ 3> |Aj| v 2 . In this 
case the heavy states H, H and A are nearly mass degenerate and decouple from the light 
state h. 



In Fig. 12 we show the allowed values at 90% CL for (TZn,TZi2,TZi 3 ) for a general scalar 
state ip® with m^o = 126 GeV, assuming that the alignment parameters qj (f = u, d, I) are real. 
We have imposed \q u \ < 2, in order to satisfy the flavour constraints for a charged Higgs below 
the TeV scale, and moreover we have set |<^,z| < 10. It is seen that the CP-odd admixture in 
the 126 GeV state has an upper bound TZ^ < 0.7, similar to that obtained in Ref. [35] within 
2HDMs of types I and II, with explicit CP violation and soft breaking of the Z2 symmetry. 



5 Summary 

The recent LHC discovery of a new neutral boson, with mass close to 126 GeV, provides for the 
first time direct information on the electroweak symmetry breaking mechanism. The current 
data is so far compatible with the SM Higgs hypothesis, although a slight excess in the diphoton 
channel has been observed by the ATLAS and CMS collaborations. This channel is particularly 
interesting since the decay of the Higgs into two photons occurs at the one-loop level and is 
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therefore sensitive to new charged particles that couple directly to the Higgs. 

As new and more precise data become available, we shall test whether the properties of 
the 126 GeV particle correspond indeed to the SM Higgs boson or they manifest evidences 
for new phenomena, perhaps signalling the existence of a much richer scalar sector. Present 
experimental errors are still large but, nevertheless, they already allow us to extract useful 
constraints on alternative scenarios of electroweak symmetry breaking. 

2HDMs constitute the simplest extension of the SM scalar sector, satisfying the electroweak 
precision tests, and give rise to interesting new phenomena through their enlarged scalar spec- 
trum containing five physical scalars. In order to avoid dangerous FCNCs, the 2HDM phe- 
nomenology has been usually particularized to a few specific implementations, based on discrete 
Z2 symmetries, which severely restrict the fermionic couplings of the scalar bosons. The most 
widely used scenario is the so-called type II 2HDM, since it corresponds to the tree-level scalar 
sector of the minimal supersymmetric SM. However, the phenomenological FCNC requirements 
can be easily satisfied imposing a much softer alignment condition on the Yukawa couplings. 
The resulting A2HDM provides a general framework to describe an extended scalar sector with 
two Higgs doublets and no FCNCs at tree level, which includes as particular cases all previ- 
ously considered 2HDM variants. It has a much larger parameter space with plenty of new 
phenomenological possibilities, such as new sources of CP violation and tunable strengths of 
the (family universal) Yukawas. Thus, it is the appropriate framework to perform an unbiased 
phenomenological analysis of the Higgs data. 

In this paper, we have analyzed the present data on the Higgs signal strengths from the 
ATLAS, CMS, CDF and D0 collaborations, within the framework of the A2HDM. Even with 
the currently large experimental uncertainties, interesting conclusions can be obtained regarding 
the preferred regions in the parameter space of the model. We have considered a variety of 
possible departures from the SM predictions, within this framework, including the effects from 
new CP-violating phases. In particular, we have searched for possible ways to enhance the 
diphoton channel while being compatible with the rest of the data. 

The measured WW, ZZ and 77 decay channels of the new boson suggest that its coupling 
to the weak vector bosons ZZ) is close to the SM one. This rules out the possibility 

of a pure CP-odd assignment for the quantum numbers of the new Higgs-like boson. A CP-even 
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scalar, either pure or with a CP-odd admixture arising from CP-violating terms in the scalar 
potential, however, can accommodate the data rather well. 

By flipping the relative sign of the top Yukawa coupling, the top-quark contribution to the 
Higgs decay amplitude into 27 interferes constructively with the dominant W ± contribution. 
This can only be realized in the A2HDM for large values of \q u \, given that g^oyy ps Qhvv ■ 
However, flavour constraints on a charged Higgs below the TeV scale (from Z — > bb, b — > 57 
and B°-B° mixing) require that \q u \ < 2, even in the most general CP-violating A2HDM. Thus, 
a 27 enhancement through a constructive interference of the top and W ± contributions could 
only be possible in a decoupling scenario with an enormously large H ± mass. 

Including the charged scalar contribution to the Higgs decay amplitude into two photons, 
one can explain the observed excess without significant deviations of the neutral scalar couplings 
from the SM limit, and satisfying at the same time the flavour constraints. This appears to 
be the most natural and likely possibility to accommodate current data within the A2HDM 
framework. The confirmation by future data of a significatively enhanced 27 decay width could 
be a strong indication that a light charged scalar is around the corner, within the LHC reach. 

The possibility that a CP-even and a CP-odd Higgs bosons have quasi-degenerate masses 
near 126 GeV was also analyzed. An excess in the 77 channel can occur in this case due to the 
contributions from both scalars (when signal strengths are added incoherently). We have also 
considered the most general A2HDM with complex Yukawa couplings. Since the Higgs signal 
strengths are CP-even observables, there is no interference between the contributions from the 
real and imaginary parts of the Yukawa couplings. It is then possible to enhance the 77 decay 
rate with a complex Yukawa coupling which has its real part close to the SM-like limit. 

Future improvements of the present bounds on neutral and charged Higgs bosons, or perhaps 
their direct discovery, as well as more precise measurements of the current Higgs signal strengths 
are expected from the LHC in the next years. The complementarity between flavour constraints 
and collider searches for new scalar resonances will be crucial for the understanding of the 
mechanism of electroweak symmetry breaking. We have shown different alternative scenarios 
within the A2HDM that can accommodate present data very well, placing bounds on the 
relevant parameter space and discussing possible consequences that could be tested in the near 
future. 
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A Scalar Potential 



In the Higgs basis, the most general scalar potential takes the form 



+ A x (V^i)' + A 2 ($ 2 $ 2 ) 2 + A 3 ($ 2 $ 2 ) + A 4 ($1$ 2 ) 



+ 



A 5 ${$2 + A 6 + A 7 $£$ 2 ) $i$2 + h.c. 



(30) 



The Hermiticity of the potential requires all parameters to be real except A 5 , A 6 and A 7 ; 
thus, there are 14 real parameters. 

The minimization conditions (0|<&f (x)|0) = (0,v) and (0|$f (x)|0) = (0,0) impose the 
relations 

Hi = -Xiv 2 , fi 3 = -^X 6 v 2 . (31) 

The potential can then be decomposed into a quadratic term plus cubic and quartic inter- 
actions 



V = --\ 1 v A + V 2 + V 3 + V i . 



(32) 



The mass terms take the form 



V 2 = M 2 H± H + H- + £ (5i,5 2 ,S 3 ) X 



with 



and 



= M|± + - /i 2 + -MfjH 2 + - M^A 2 

2 2 2 



M 2 ± = /x 2 + - A 3 u" 



(33) 



(34) 



/ OX ..2 



A< = 



-v 2 \\ 



\ 



-v 2 \\ 



(35) 



2A^ 

v 2 \f M 2 H± +v 2 (^ + \f) 

-v 2 \\ -v 2 \\ M 2 H± + v 2 (4f - A*) ; 

where Af = Re(Aj) and \] = Im(Aj). The symmetric mass matrix .M is diagonalized by an 



V 
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orthogonal matrix 1Z, which defines the neutral mass eigenstates: 

( Ml \ / h \ 





m = n 1 



\ 



Mjj 



Ml j 



H 



R 



fs x \ 

s 2 

\S 3 J 



(36) 



Since the trace remains invariant, the masses satisfy the relation 

Mi + M 2 H + M\ = 2M 2 H± + v 2 (2Ai + A 4 ) 



(37) 



The minimization conditions allow us to trade the parameters /xi and /J3 by v and A6- The 
freedom to rephase the field $2 implies, moreover, that only the relative phases among A5, \$ 
and A7 are physical; but only two of them are independent. Therefore, we can fully characterize 
the potential with 11 parameters: v, /12, (Ai,...^!, arg(AsAg) and arg(A5A^). Four parameters can 
be determined through the physical scalar masses. 

In the CP conserving limit A5 = Ag = A^ = and £3 does not mix with the other neutral 
fields. The scalar spectrum contains then a CP-odd field A = S3 and two CP-even scalars h 
and H which mix through the rotation matrix pi). In this case, the scalar masses are given by 



Ml 



where 



A) 



Ml 



:s+a) , 



Ml 



M 2 H± + 



A4 
2 



A 



R 



M 



2Ai 



A 



Ml ± + v 2 



-2X l + ^ + \f 



and the mixing angle is determined through 



tan a 



M 2 - 2\ lV < 



(38) 

(39) 
(40) 

(41) 



We use the notation M w o to emphasize that these are the neutral scalar masses in the CP- 
conserving limit. The cubic and quartic Higgs couplings involving the charged and the neutral 
physical scalars (without Goldstone boson couplings) take the form, 

V, = vH+H- (x 3 s 1 + xfs 2 -x I 7 s 3 )-lvx I 7 si-lvx I 7 sis 3 -lvx I 6 sls 3 



+ X 1 vS 3 1 + ^vX 7 l Sl + ^vXfS 2 1 S 2 + ^v (2A^ + A 3 + A 4 ) S 1 S [ 
-\v (2A^ - A 3 - A 4 ) S X S 2 + \vX\ x S 2 S 2 3 -2vX{ S 1 S 2 S 3 , 



(42) 
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Vi = H + H (^)\2H + H + S 2 + A2 5*1 + A2 5*2 — Ay S1S3 + A R SiS^j 
+ \ (A 3 + A 4 + 2A*) (S^) 2 + l - (A 3 + A 4 - 2A*) (S^Y + ^ (S 2 S 3 



1 x 1 

_ \I C3c \Jc2c C 7 

2 ' 
Ai 
4 



2 



A R 

\1 o3 G \J C C '"To c2 C c o3 i 6 o3 G i /y 7 C c3 i /v 7 rf o o2 
Ag O-LO3 — A 5 1 02»J3 — Oi0 2 >J3 — 0\0 3 H — 1 02 H — OiD 2 i — O1O2O3 



A? 
2 



A? 
2 



(43) 



In the CP-conserving limit all vertices involving an odd number of 5*3 fields vanish. 



A.l Neutral scalar mass matrix to lowest order in CP violation 



Assuming that A5 and A 6 are small, we can diagonalize the mass matrix (35) perturbatively 
as an expansion in powers of these CP-violating parameters. The leading corrections to the 
neutral scalar masses are quadratic in A^: 

M 2 o = M 2 + af (A5) 2 + af (Ag) 2 + of , (44) 



where M o denote the corresponding masses in the CP-conserving limit given in (38) and 



a: 1 



(MJ - 2A lW 2 ) 



2\ lV 2 + M 2 -M 2 H - Mi 
2v 6 \* 



(45) 



The physical states tp® = {h, H, A} receive corrections at first order in Al, 6 , which are given 



by 



/ h \ 
H 



^ cos a sin a ei 3 \ f 
— sin a cos a e 2 3 S2 
£31 £32 1 J \ S 3 ) 



(46) 



where 
ei3 = 

£31 = 



(sin a A5 + cos a A 6 ) 



{M\ - M 2 ) 
"2^2 KA I 5 + 2a 2 4 A I 6 ) 



£23 



^32 



{M\ - M 2 H ) 
I 



(cos a \\ — sin a A 6 ) 



2v 2 



(47) 
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Note that for the case of a scalar potential with a softly-broken Z 2 symmetry in the Higgs basis 
we have A 6 = A 7 = and, therefore, e 3i = 0. 



B Scalar Couplings to the Gauge Bosons 

The scalar doublets couple to the gauge bosons through the covariant derivative and gauge- 
fixing terms: 

2 

C K + ^2(D^ a y D^ a + C GF = C V 2 + + C^v + C<j? v + + Cpv* , (48) 

i=l 

where Ck is the usual gauge-boson kinetic term and the covariant derivative is given \rj^\ 

Dn = d u + ieQAn + i — — — ZJT 3 — Q sin 2 9w) + W \T+W}, + TL Wj . It is convenient to adopt 

COS W >~ 

the following gauge-fixing term (£ = 1), 



C 



(;F = -- (d^f - - (d^ + M Z G°) 2 - (d»Wl + iM w G + ) {d v W - iM w G-) , (49) 



which cancels exactly the quadratic mixing terms between the gauge and Goldstone bosons 
generated by the covariant derivatives, so that C^y = 0, and provides the Goldstone bosons 
with the masses M G ± = My/ — gv/2 and M G o = Mz = Mw/ cos^vy. Then, 

£v* = -\ (<V) 2 - \ ^ z "f + \ M l z ^ - ( d "K) ( d » wv ) + M w K wlM > ( 5 °) 

while 



+ \ d^G° d^G - l - M\ (G ) 2 + <9 M G+ d^G" - M w G + G~ 
The interaction terms between the scalar and gauge bosons are given by: 
L^y = ie [A^ 1 + cot (29 w ) Z M ] \(H + d^H~) + (G + d^G~ 



+ 



sin (26 



w , 



(G°d,S 1 ) + (S 3 d»S 2 ) 



+ 9 -W^ [(H-d l ^S 3 )-z(H-d,S 2 ) + (G-d^G°)-z(G-d,S 1 
+ | W [(H + d, S 3 ) + 1 (H + d, S 2 ) + (G+ d» G°) + i (G+ d, Si) 



(51) 



(52) 



^The weak mixing angle 9w is denned through the relation gsmQw = g' cos 0w = e. The operators T± = 
4=(Ti ± T-i) and T3 can be expressed in terms of the Pauli matrices by Tj = , 
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c 



4>V 2 



Si 



-MlZ^ + M^WlW* 



+ (eM w # -gM z sm 2 e w Z») (G+W^ + G~W§ 



(53) 



£ 



4P-V 2 



+ 



-M 2 z Z^ + M^WlW» 



[H 2 + h 2 + A 2 + {G Q f] 



+ je 2 [A" + cot (26 w ) Z^f + ^ WlwX (G + G- + F + iT) 



{A<* - tan^ ^ (G+W^ + G~Wl) + S 2 (H+W^ + H~W§ 

+ z5 3 (^-^t-if+^)+ ? G' (G'-^t-G' + ^)l, (54) 



with Si = Kjitf ((p°j = {h, H, A}) and the usual notation A d M B = A(d^B) - (d^A)B. 



C Statistical treatment and data 



To obtain the preferred values for the parameters of the A2HDM we build a global \ 2 function 



X 



(/j>k - faY 



at 



(55) 



where <7j is the experimental error extracted from the data at 1 a. Errors on the reported Higgs 
signal strengths fa are symmetrized using 



8 fa 



{Sfi + ) 2 + (8 fa 



(56) 



where 8jl± are the one-sided errors given by the experimental collaborations. We use the latest 
data available after the "Hadron Collider Physics Symposium 2012 (HCP2012)", including 
the latest update from ATLAS of the high-resolution channels 77, ZZ^*> 111. For the diphoton 
channels we use the data given by ATLAS and CMS at 7 and 8 TeV, provided in Refs. 1 2 p3p4 



For the rest of the channels we use the averages listed in Table |2j which include the 7 © 8 TeV 



data reported by ATLAS and CMS together with CDF and D0 data 10-12 at s/s = 1.96 TeV. 



For a general channel with inclusive production we have (neglecting the subdominant pro- 
duction channels) 

(57) 



°gg Br(y° -> fe) 
off ' Br(^ fc) SM ' 
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Table 2: Higgs signal strengths in each of the channels considered in this work. Averages obtained 
from ATLAS and CMS data at 7 8 TeV together with CDF and D0 data at y/s = 1.96 TeV. (*) 
We do not consider non-inclusive measurements in the tt channel. Due to the large current errors 
associated with these measurements, our conclusions would not be modified at this level. 



Channel 


h 


Comment 


bbV 


1.1 ±0.44 


ATLAS, CMS, CDF and D0 


1,2 


11 


12 


( 


'our average) 


WW]] 


-0.2 ± 1.56 


ATLAS and CMS |l| 


y 


14 


(our 


average) 


WW 


0.76 ±0.21 


ATLAS, CMS, CDF and D0 


III 


2,10 


11 


13 


(our average) 


zz 


0.96 ±0.26 


ATLAS and CMS 


15 


I (our 


average) 


tt (incl.) (*) 


0.89 ±0.86 


ATLAS and CMS 


15 


I (our 


average) 


77 


1.66 ±0.32 


ATLAS and CMS 


15 


I (our 


average) 


1133 


2.18 ±0.84 


ATLAS and CMS 


15 


I (our 


average) 



For the Higgs searches in the 77 channel, the ATLAS and CMS collaborations have established 
different categories. To take this into account, we write the Higgs signal strength in a given 77 
channel as 

<p1 tggF CTggF + CyBF VyBF + ^VH OyH Br(y?° -> 77) 

e ggF VggF + e VBF <7 V BF + e VH O'vH ±jr lV 9 i ~* 77JSM 

where the coefficients e( 95 F,vBF,VH) accounting for the relative weight of each production channel 
have been provided by ATLAS and CMS (5l|6j. The top-quark-fusion contribution could be 



added in a similar way. In Eq. (58), the SM production cross sections and decay widths are 
taken from the web page of the LHC Higgs Cross Section Working Group [80]. For the gluon- 
fusion production mechanism we have 

a(gg^tf) = a ggF = cgo*f F , (59) 

<p° n 

where the scaling of the gluon-fusion cross section C g g was defined in section 3 Vector-boson 
fusion scales with the coefficient Tin as 

a(qq' ^qq'rf) = a YBF = (TZ a ) 2 a^ F , (60) 

and similarly for the associated production with a vector boson 

*(qq^Vtf) = a VH = (^i) 2 4h- (61) 
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Figure 13: Diagram contributing to the one-loop ip®H + H vertex correction. 

D Perturbativity Constraints 

The charged Higgs boson contribution to (p® — > 77 depends crucially on the value of the neutral 
scalar coupling to a pair of charged Higgs bosons. To assure the validity of perturbation theory, 
upper bounds on the quartic Higgs self-couplings are usually imposed requiring these to be 



smaller than 8ir (see 30 31 and references therein). The cubic Higgs self-couplings are also 
bounded indirectly in this way. In this work we consider an alternative perturbativity bound 
on the relevant Higgs cubic coupling which is more restrictive for light charged Higgs masses. 



Consider the f>^H + H~ one-loop vertex correction given by Fig. 13 The contribution of this 
diagram is finite and can give us an idea about the allowed magnitude of the cubic coupling in 
order not to spoil the perturbative convergence. We obtain: 



°H+H~)eS 



where 



7> 2 A 2 

V \°H+H- 

16tt 2 M^ ± 



1-2/ 



Z 



M 2 H± 



\9H+H~ 



;i+a) 



(62) 



Z{X) = I dy I dz [(y + z) 2 + X(l-y-z-yz)] \ (63) 
Jo Jo 

Allowing the correction to be at most 50% (A ^ 0.5) constraints the allowed parameter space 
in the (\^o H+H - , M H ±) plane to be within the blue (hashed) region indicated in Fig. |6| 



E Oblique Parameters 

Possible deviations from the SM in the gauge-boson self-energies are usually characterized 
through the oblique parameters S, T and U 81 . Taking as a reference SM Higgs mass M^ re f = 



126 GeV, the most recent global fit to electroweak precision observables quotes the values 67 
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68 : 



S = 0.03 ±0.10, T = 0.05 ±0.12, U = 0.03 ± 0.10 . (64) 

The expressions for the oblique parameters in the CP conserving A2HDM are adapted from 



Ref. 82 : 



S 



ttM 2 



cos 2 a 



+ sin 2 a 



B 22 (M 2 Z ; Ml M 2 h ) - M 2 Z B (M 2 Z ; M|, M 2 ) + B 22 (M 2 Z ; M 2 H , M 2 
B 22 (M 2 Z ; Mf, M%) - M 2 Z B (M 2 Z ; M|, M 2 h ) + B 22 (M 2 Z ; M 2 , M 2 



B 22 (M 2 ; M 2 H± ,M 2 H± )- B 22 (M Z ; M z , M 2 ) + M 2 B (M 2 ; M 2 , M 2 ) } , (65) 



T 



cos 2 a 



+ sin 2 a 



F(M 2 H± ,M 2 H ) - T{M%, M 2 ) + 3F(M 2 Z , M 2 ) - 3^(M^, M 2 ) 
7-(M| ± , M 2 ) - 7 {Ml M\ ) + 3 T{M\, M 2 H ) - 3 T{M 2 W , M 2 H 



+ F{M 2 H± ,M 2 A ) - 3F(M Z , M 2 ) + 3 T{M^, M 2 



(66) 



U = H{Ml)-H{M 2 z ) + 



where 



TlM 2 , 



sin 2 a B 22 (M^; M 2 H± , M 2 ) + cos 2 a B 22 (M^; M 2 H± , M 2 



+ B 22 (M^;M 2 H± ,M 2 A )-2B 22 (M^M 2 H± ,M 2 H± ] 



TlM 2 



n(M 2 ) 



TlM 2 



cos 2 a 



+ sin 2 a 



sin 2 a B 22 (M 2 ; M 2 , Mj) + cos 2 a B 22 (M 2 ; M 2 H , M 2 A 
- B 22 (M 2 Z ;M 2 H± ,M 2 H± ] 

B 22 (M 2 ; My, M 2 h ) - M v B (M 2 ; My, M 2 
B 22 (M 2 ; My, M 2 H ) - My B (M 2 ; M 2 , M 2 H 



(67) 



B 22 {M 2 - M 2 , M 2 ) + M v B (M 2 ; M 2 , M t 



h,ref) 



(68) 
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The loop functions are given by 



B 22 (q 2 ; m\,m\) = i (A + 1) \m\ + ™l ~ \ Q 2 ] ~ \ J dx X log (X - it) 

B (q 2 ;ml,ml) = A— / dx log (X — ie) , 

Jo 



(69) 
(70) 



2 2 A 1 / 2 . 2x m\m\ ( ml 

' - -(m 1 +m 2 \) 2 log — 

2 mf — m 2 \ m 2 



with 



F(jni,mi) = o K + mi) - _ 2 ^_ 2 log ( ^ ) , (71) 

2 

X = m\x + ml(l -x) -q 2 x(l - x) , A = - — - + ln47r - j E , (72) 



in d space-time dimensions, where 7^ is the Euler-Mascheroni constant, and where we have 
defined: 

B 22 (q 2 ;ml,m 2 ) = B 22 (q 2 ; ml, ml) - B 22 (0; ml, ml) , (73) 
B {q 2 -mlm 2 2 ) = B (q 2 ; m\, m 2 2 ) - B (0; m\, m 2 2 ) . (74) 
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